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MSSM is excluded. The minimal supersymmetric fermiophobic Higgs scenario can naturally 
be formulated in the context of the NMSSM that admits Z 3 discrete symmetries. In the 
O fermiophobic NMSSM, the SUSY naturalness criteria are relaxed by a factor N c yf/g A ~ 25, 
removing the little hierarchy problem and allowing sparticle masses to be naturally of order 2- 
3 TeV. This scale motivates wino or higgsino dark matter. The SUSY flavour and CP problems 
as well as the constraints on sparticle and Higgs boson masses from b — > 37, B s — > /i/i and direct 
LHC searches are relaxed in fermiophobic NMSSM. The price to pay is that a new, yet unknown, 
mechanism must be introduced to generate fermion masses. We show that in the fermiophobic 
NMSSM the radiative Higgs boson branchings to 77, 7Z can be modified compared to the 
fermiophobic and ordinary standard model predictions, and fit present collider data better. 
Suppression of dark matter scattering off nuclei explains the absence of signal in XENON100. 



ABSTRACT 



If a light Higgs boson with mass 125 GeV is fermiophobic, or partially fermiophobic, then the 
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1 Introduction and motivation 



The TeVatron [JJ |2] and LHC experiments presented their new and updated results [3] IU G2, EJ 
El El EH El Q2] on searches for the Higgs boson [B] at the Moriond 2012 conference [H]. 
While on average the data is consistent with the standard model (SM) Higgs boson with mass 
125 GeV, interesting anomalies start to emerge that may signal unexpected new physics in 
the Higgs sector. The most interesting of them is a local 3a level excess in searches for the 
fermiophobic (FP) Higgs boson [15] in 77 final states both in the ATLAS and CMS experiments 
[16J. This signals that there is an anomalously large contribution in the observed 77 excess 
coming from the vector-boson fusion (VFB) Higgs production mechanism. Indeed, the relative 
weight of the latter and the associate production with W, Z (VH) is enhanced with respect to 
the SM dominant gluon-gluon fusion channel (ggF) in the FP high-p^ selections applied by the 
CMS and ATLAS. This anomaly is accompanied by a deficit of WW* compared with the SM 
in all experiments. 

The Higgs boson mass w 125 GeV is peculiar in several ways. In the context of FP Higgs 
boson, there is an accident that at the 7-8 TeV LHC the Higgs boson signal rate in the 77 
channel, a x BR, happens to be equal to the SM one in the vicinity of this Higgs mass value [T7j. 
This is the reason why the LHC inclusive 77 excess is consistent with the FP Higgs boson. At 
the same time, the signal rates for other gauge boson channels, WW*, ZZ* , Zj, are predicted 
for be 40-50% suppressed compared to the SM [TTJ . The M h m 125 GeV Higgs boson is peculiar 
also because, for a SM-like Higgs of that mass, branching fractions for many decay channels 
are measurable at the LHC. Therefore, the LHC is, in principle, able to determine the nature 
and properties of the 125 GeV Higgs boson. 

Motivated by these results we performed a global fit to all available collider data to determine 
which Higgs boson scenario is currently favoured [18] . improving and extending similar pre- 
Moriond fits [19] . A purely FP Higgs boson gives a fit to present data almost as good as the 
SM one, but with very different predictions for the signal rates at the LHC. A partially FP 
Higgs boson, however, gives a significantly better fit to current data than the SM [T5]. This 
is because the FP Higgs qualitatively describes the observed anomalies in the data correctly, 
although it predicts larger signal rates than observed. Small additional branching fractions into 
fermionic channels, for example into bb, decrease the overall rate and improve the fit to data 
significantly. 

Partial fermiophobia is exactly what is expected to happen when considering the FP Higgs 
boson scenario as an effective low energy theory in the context of quantum field theory [20J. 
Because the SM fermions are massive, at loop level non-vanishing Yukawa couplings are gener- 
ated even if at some high new physics scale the Higgs boson was initially purely fermiophobic. 
As long as the fermion mass generating mechanism is unknown, the induced Yukawa couplings 
represent an uncertainty in the low energy FP Higgs model. Present data suggest that the 
Higgs boson might be partially fermiophobic [TTJ |20l [21] . 

The 125 GeV Higgs boson is also unsatisfactory, because this mass is below the SM vacuum 
stability bound in case new physics appears only above the scale of gauge coupling unifica- 
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tion [22]. The vacuum stability can be made consistent with the 125 GeV Higgs mass [23] by 
extending the scalar sector with dark matter candidates [2H [25] . But the simplest solution is 
provided by (partial) fermophobia. Because the vacuum instability is caused by the large top 
Yukawa coupling in the SM, reducing its value makes the hinted Higgs mass compatible with 
the GUT scale. 

It is well known that the existence of the Higgs boson rises a question of why the electroweak 
scale is so much smaller than the Planck scale. The most elegant solution to that problem is 
given by supersymmetry (SUSY). However, direct and indirect collider bounds, cosmological 
dark matter abundance and constraints from dark matter direct detection experiments together 
with a Higgs boson mass of 125 GeV impose stringent constraints on SUSY scale in most popular 
SUSY models - the MSSM [261 [27] and NMSSM [281 [29]. The constraint M SU sy > 1 TeV > M z 
reintroduces severe fine tuning to theory, known as the little hierarchy problem [30], that makes 
SUSY as a solution to the hierarchy problem unnatural. 

If the Higgs boson turns out to be fermiophobic, some SUSY models are in even more serious 
trouble. We show that fermiophobia and a Higgs boson mass at 125 GeV together exclude all 
versions of the MSSM independently of any model detail. This is because in the MSSM the 
upper bound on the tree level Higgs boson mass is Mz, and large radiative corrections, domi- 
nated by stop contributions, are needed to reach 125 GeV. Fermiophobia removes the dominant 
stop loops as they are induced by Yukawa couplings. There might be large trilinear scalar cou- 
plings, the A-terms, but their contribution to M\ is negative. Dimensionful trilinear couplings 
may trigger electroweak symmetry breaking [3T] via a dimensionful Coleman- Weinberg [32J 
mechanism, but cannot increase the Higgs boson mass. 

Nevertheless, SUSY models with additional tree level contributions to the Higgs boson mass, 
such as the NMSSM, are viable fermiophobic SUSY theory candidates. In fact fermiophobia 
can cure some SUSY problems and make it more compatible with experimental data. Firstly, 
the fine tuning problem of SUSY, also coming from loop contributions to the Higgs boson 
mass squared, is now induced by gauge couplings, improving the fine tuning by a factor of 
N c yf/g A ~ 25. This improvement completely removes the little hierarchy problem. SUSY 
masses of order 2-3 TeV become completely natural. Allowing for some fine tuning, even the 
split SUSY [33] with a very heavy scalar sector becomes viable. This would explain why no 
sparticles have been discovered by the LHC up to now [34, 35J. Secondly, SUSY flavour and 
CP problems are improved by removing (or decreasing) the Yukawa couplings and by allowing 
also squark and slepton masses to be at a few- TeV scale. Thirdly, the additional Higgs bosons 
can be light. For example, the absence of constraints from 6—7-57 an d B s — > \x\i [30] allows 
the charged Higgs boson to be light, opening again the possibility for its discovery at the LHC. 
Fourthly, the SUSY fermion sector may be either light or heavy. While gluinos could have been 
abundantly produced at the LHC if they were very light, for colourless fermions there would 
exist only collider bounds from LEP and Tevatron. Fifthly, the constraints on dark matter 
are relaxed. Because neutralino elastic scattering off nuclei is dominated by tree level Higgs 
boson exchange, this process would be suppressed and the prospects for dark matter discovery 
at the XENON100 are decreased in this scenario [37]. As the SUSY scale could now be large, 
higgsino or wino relic abundance would become a natural explanation to the dark matter of 
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the Universe. 



The obvious question in any FP Higgs boson scenario is what is the alternative mechanism for 
generating the observed fermion masses. Because the top quark mass is so large, it cannot be 
generated radiatively. The most plausible scenario for generating such large fermion masses 
is strong dynamics above the electroweak scale [38J. In such a scenario both the composite 
Higgs boson fermiophobia and fermion masses might originate from the same new physics. 
A generic prediction of strong electroweak symmetry breaking scenarios, including composite 
Higgs models, is the appearance of new resonances at 2-3 TeV. In the following, we assume 
that such or any other new physics scenario above the electroweak scale generates the top quark 
mass. 

In this work we formulate a FP NMSSM as a minimal FP SUSY model. Originally the NMSSM 
was constructed to solve the fi problem of the MSSM and to have an additional contribution to 
the masses of the Higgs bosons. To achieve that, an additional Z 3 symmetry is usually imposed 
on the NMSSM. We show that choosing quantum numbers of this symmetry appropriately, 
superpotential Yukawa terms may be forbidden in the NMSSM. In this case, all the pros and 
contras of FP SUSY discussed above apply also to the FP NMSSM. In our phenomenological 
study of the model we concentrate on radiatively induced Higgs boson decays h —> 77 and 
h — > Z7. This choice is motivated by the fact that SUSY effects most easily show up in loop 
level processes. Our aim is to study whether the FP SUSY Higgs boson gives a better or worse 
phenomenological fit to the LHC data than the FP SM Higgs boson. We find that the new 
SUSY contribution can enhance or reduce the 77 and Zj rates as much as 50% compared to 
the FP SM Higgs depending on the sign of the fi parameter. Because the data prefers smaller 
rates [IS] , the FP SUSY Higgs can give a better fit to data than the SM. 

We note that in the MSSM the WW* and ZZ* rates can be reduced at tree level compared to 
the SM. At the same time, the decay h — > 77 is dominated by the VT-boson loop, introducing 
a correlation between the two processes. Should the observed deficit in WW* persist together 
with the 77 excess, new physics beyond the MSSM would be required. A FP NMSSM might 
be a good candidate for such new physics. 

The paper is organized as follows. In section 2 we study fermiophobic Higgs scenarios in SUSY. 
In section 3 we study the radiative decays h — > 77 and h — > Z^y in a specific parameter region 
of the FP NMSSM model. In section 4 we add some discussion and conclude in section 5. 



2 Fermiophobic supersymmetry 

To formulate a supersymmetric FP Higgs boson theory, the first attempt should be made in the 
MSSM. However, we are going to show that the Higgs boson mass w 125 GeV is by far too 
large to be generated in the FP MSSM since loop corrections from the top Yukawa coupling 
are absent. The dominant SUSY loop contribution to the Higgs mass in the FP MSSM comes 
from the large trilinear A-term, but this contribution is always negative. Thus the FP MSSM is 
definitely excluded on phenomenological grounds. In order to rescue supersymmetry, we show 
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that the NMSSM offers a natural framework to formulate a supersymmetric FP Higgs scenario 
consistent with experimental results. 



2.1 Fermiophobic MSSM 

The well known MSSM superpotential is 

W = y u QH u u c + y d H d Qd c + y e H d Le c + fiH u H d , (1) 

where y u , y d and y e are the up quark, down quark and charged lepton Yukawa couplings. The 
LHC hints for a FP Higgs imply that at least the third generation Yukawa couplings must be 
strongly suppressed compared to their SM values so that the production mechanism gg — > h 
and the decay channels h — >■ bb, h — >■ rf become subdominant compared to the gauge boson 
processes. Following this indication, in this paper we study its implications without trying to 
explain the suppression of Yukawa couplings in the context of MSSM. Thus, for simplicity, we 
just take y u = y d = y e = 0. 

In the MSSM the tree level Higgs boson mass has the well known upper bound M^ ree < Mz- 
This comes from the fact that in SUSY the Higgs quartic couplings are generated by gauge 
couplings via the D-terms. As the Higgs boson quartic coupling is the only free parameter in 
the SM Higgs sector, in the MSSM there is no freedom to tune the tree level Higgs boson mass. 
To be consistent with experimental data, in the MSSM very large positive loop corrections to 
M\ must be generated. Those loop corrections are dominated by top squark contributions that 
are induced by the top Yukawa coupling y t 



AM 2 = 3y 4 



v sin /3 
8tt 2 



Ml , XI / X 2 



l oe; — 5. _i_ , 
8 mf Mf V 12M b .i 



(2) 



where Ms is the average stop mass, tan (5 = v u /v d , v 2 = v^ + v^, and X t is the stop mass mixing 
parameter 

X t = A t — fi cot — — — a cot 0, (3) 
Vt 

where a t is the trilinear coupling of the soft term a t QH u u c . In order to achieve the Higgs boson 
mass indicated by the LHC experiments, rj 125 GeV, Eq. (|2]) implies that the stop masses 
must exceed TeV scale. 

In the FP MSSM the dominant stop contribution is absent since we take y t — > 0. However, 
loops induced by very large trilinear soft interaction a t QH u u c with a t > 1 TeV are still allowed 
in general MSSM even in the absence of Yukawa couplings. In the FP limit we obtain 

. .^n 3v 2 sin 4 /3 af . 

AM, = — — - — (4) 

h 8vr 2 12M|' 1 ] 

that is negative. The dimensionful couplings like the A-terms may be used to generate negative 
Higgs mass terms radiatively, thus generating dynamical breaking of electroweak symmetry [31 J, 
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but they do not increase the Higgs mass prediction in the MSSM. This is because the Higgs 
boson quartic self coupling is fixed by the D-term that gives the upper bound. As the chargino 
loop contributions to the Higgs boson mass are of order few GeV, in the FP MSSM the Higgs 
boson mass 125 GeV is not achievable. Independently of model details, the FP MSSM is 
excluded by the Higgs boson mass. 



2.2 Fermiophobic NMSSM 

NMSSM is the next to minimal supersymmetric standard model whose particle content is 
extended by a gauge singlet chiral superfield S (for reviews and references therein see [22] )• The 
original motivation for the NMSSM was to explain why the MSSM superpotential parameter 
fiH u H d is of the same order as the soft SUSY breaking parameters. In addition, in the NMSSM 
the Higgs bosons obtain tree level mass not determined by the D-terms, thus allowing larger 
Higgs masses than Mz- To achieve those goals, usually the most general NMSSM is constrained 
by imposing an additional Z% symmetry in addition to the i?-parity. Those properties make 
the NMSSM our prime candidate for the minimal FP SUSY model. The superpotential of the 
FP NMSSM that we would like to obtain is given by 

W = \SH u H d + ^S 3 , (5) 

together with the following soft SUSY breaking terms 

£ so ft = - {m\ji\h u + m\h ] d h d + m^sj - (^a x sh u h d + ^a k s 3 + h.c.J , (6) 

where s stands for the scalar component of the singlet chiral superfield S. Thus we have to 
forbid the fj,H u H d , S, S 2 , y u QH u u c , y d H d Qd c , y e H d Le c terms by imposing an additional Zn 
symmetry and appropriately choosing the charges to satisfy the following constraints 

(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

One could choose X L = Xq and X e c = X d c, then ^ would be superfluous. 

For Z 3 , if Xh u = Xh a = Xs = 1 and Xf ermion = 0, the above equations are satisfied. In addition, 

the lowest order Yukawa couplings can be generated via d = 6 operator d ^ 2 d QH d d c , where A 
is the scale of new physics. This demonstrates that small but non-vanishing Yukawa couplings 





+ x Hu 


+ X u c 


7^0 


mod N, 




> + X H d 


+ X d c 


7^0 


mod N, 


x L + x Ha 


, + X e c 


7^0 


mod N, 




Xh u 


+ x Hd 


7^0 


mod N, 






X s 


7^0 


mod N, 






2X S 


7^0 


mod N, 






SX S 


= 


mod N, 


X s 


+ x Hu 


+ x Hd 


= 


mod N. 



5 



should exist also in the FP Higgs scenarios. On the other hand, if Xjj n = Xjj d = X$ = 1 and 
-^fermion = 2, the above equations are also satisfied, but the lowest order Yukawas would be 
generated as e.g. ^-QH d d c . These charges satisfy two additional equations 



This could possibly generate the Yukawa couplings for the first two generations but not for 
the top quark. Therefore we have to assume that the significant amount of third generation 
fermion masses should come from some additional mechanism. The prime candidate for such 
a mechanism is some strong dynamics above 2-3 TeV scale. 

The Z 3 symmetry could come from the breaking of an U(l)' [40J. In this case one has to 
keep in mind the possibility of discrete gauge anomalies |41j . In the NMSSM these have been 
considered in [12]. If one chooses -Xf ermion = 0, the anomaly constraints can be simply evaded. 

Breaking of the Z 3 symmetry in the early universe could create a problem with the domain 
walls that produce an anisotropy in the CMB and ruin nucleosynthesis [43J. The problem can 
be solved by allowing for radiative generation of small renormalisable Z 3 breaking terms jH] . 



3 Phenomenology of FP NMSSM Higgs bosons at the 



Our aim in this section is to study radiatively induced decays h — > 77 and h — > Z*y of the 
SM-like FP Higgs boson in the context of NMSSM. If there are light superparticles or light 
additional Higgs bosons, their effects are first expected to show up in loop level processes. 
However, before proceeding with this study, we have to show that the FP NMSSM is a viable 
model. Therefore we start studying the FP NMSSM Higgs sector. We do not attempt to scan 
the full parameter space of the model. Instead, we start by fixing the parameters in the Higgs 
sector to one particularly interesting point with tan/? = 1 and decoupled CP-even singlet, that 
allows for two following Higgs scenarios: (i) the 125 GeV excess is due to the lightest CP-even 
Higgs boson, the remaining neutral Higgs bosons are too heavy and have no direct decays to 
two gauge bosons; (ii) the 125 GeV excess is due to the next-to-lightest CP-even Higgs boson, 
the lightest one has no direct decays into two gauge bosons and remains invisible at the LHC. 
Here, we will focus on scenario (i). Then, we relax the condition tan/3 = 1 and analyze the 
impact of a tan (3^1 value, in the approximation in which the scalar singlet is very heavy. 
We finally show how the FP NMSSM can be distinguished from the FP SM by studying the 
radiatively induced Higgs boson decays. 



Xq + X Hu + X u c + X s — 0, 
X Q + X Hd + X d c + X s = 0. 



(15) 
(16) 



LHC 
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3.1 Scalar potential and masses 

The FP NMSSM scalar potential derived from Eq. (|5]), Eq. ^ and also from the D-term 
contributions is 



V 



« + |A*| 2 ) (|^| 2 + |^| 2 ) + « + |A.| 2 ) (\hf + |^| 2 ) 
+m 2 s \s\ 2 + (a x (KK ~ h °u h d) s + ^s 3 + h.c.) 



+ \\{Kh d -h°X) + k 



s 2 \ 2 



+^ (Kf + \K\ 2 - \hf - \h-\ 2 ) 2 + 1 \h+h° d * + hlh-*\ 2 • (17) 



We suppose for simplicity that all the parameters in Eq. (17) are real. We have checked that 
the potential is always bounded from below. Supposing that only the real parts of Higgs bosons 
can get vacuum expectation values (VEVs) different from zero and parametrizing the fields as 

H ° d = 71 ^ d + H ° dR + lH ° dl) ' h ° u = 71 ^ u + H ° uR + lk ° uI ^ ' s = ^{v s + s R + is Z ) , 

(18) 

we get the following equations for the stationary points 

- 4v s sin (V2a x + kX 2 v s ^j + cos /3 (-v 2 sin 2 (5 (g{ + g\ - 4A 2 ) + 8m 2 hd + AX 2 v 2 s ) 
Wco^(3{gl + gl)^ =0, (19) 

- Av s cos (y2a x + kX 2 v s ^j + sin (3 (v 2 sin 2 (3 (g\ + g\) + 8m 2 hu + AX 2 v 2 s ) 
-w 2 sin /3 cos 2 /3 (g\ + g 2 2 -AX 2 ) ) = 0, (20) 

v s (V2a k v s + X 2 (2k 2 v 2 s + v 2 ) + 2m|) - -v 2 sin(2/3) [V2a x + 2kX 2 v s ^j = 0. (21) 

First of all we must avoid that v u — — v$ — is a minimum of the potential. This can be 
done by requiring 

m 2 h m\m 2 s < 0. (22) 

Now let us give a look at the Hessian matrix in the minimum. In the basis (h dR , h° aR) sr) we 
have for the CP-even Higgs bosons 

M| n M| 12 Mj 13 

m 2 =\ ..: mi; 22 mi; 23 | , ( 23) 
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where 
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(24) 
(25) 
(26) 
(27) 
(28) 
(29) 



So far the results have been general. However, we can see from the mass matrix that there is 
a choice of the parameters that allows no mixing between s and d , 



tan j3 


= 1, 


(30) 


k 


= A, 


(31) 




= 0. 


(32) 



Notice that tan = 1 is allowed in this model because no constraints occur from the scalar 
potential minimization nor from the Yukawa sector. Therefore this choice is the most natural 
one. Adopting, for simplicity, the choice in eqs. (30)-(32) and requiring, of course, v u ^ 0,v d 
0,vs 7^ 0, the minimization equations read 

(4m 2 d + AV) = 0, (33) 
(4t< + AV) = 0, (34) 

a ^ + m 2 s v s + X 2 vl = 0, (35) 

leading to the CP-even neutral Higgs boson mass matrix 

/ \{M 2 +v 2 \ 2 ) \{{v-vs){v + vs)\ 2 -M 2 ) \ 

M 2 = |((«- v s ){v + v 8 )\*- Ml) \ (Mf + v 2 s \ 2 ) ,(36) 

\ 2v 2 s X 2 + ^ J 

where M§ = | (g\ + g 2 ) v 2 . The corresponding eigenvalues 

A 2 ?; 2 

Ml* = , (37) 

M 2 o = M 2 + 1 -X 2 (2v 2 -v 2 ), (38) 
M 2 = ^§ + 2A 2 ^, (39) 
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correspond to the eigenvectors 



hi 
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71 
i 

71 



{Kr + Kr) , 

(h°dR ~ ^ur) ■ 



(40) 
(41) 



Two distinct phenomenologically viable scenarios occur. Since we would like to identify one 
of the CP-even eigenstates with the 125 GeV resonance hinted at by the LHC, it has to be 
made of doublets. If M 2 < M 2 then, in the notation of the MSSM, h = h\, H = hi, and 
the Higgs mixing angle is given by a = — 7r/4 = /3 — tt/2. In that case H does not have any 
direct tree level coupling to WW and ZZ that explains why the LHC does not see presently 
any other resonance, but the lightest one at 125 GeV. On the other hand, if M 2 o > M 2 o then 

h — hi, H — h^, and the Higgs mixing angle isa = 7r/4 = /3. In this case the LHC observed 
the second heaviest CP-even state because the couplings of the lightest one to fermions and to 
gauge bosons are strongly suppressed. Discovering such a light "sterile" Higgs boson is very 
difficult at the LHC. 



The CP-odd Higgs boson mass matrix in the basis (h dI , h° uI , si) is given by 
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,|A 2 
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vvs A 2 


V2 


V2 



wgA 

y 2 X 2 _ 3^ 
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(42) 



and the corresponding eigenvalues are 



M. 



Q0 



2A 2 (v 2 + v 2 s ) - 3V2a k v s - y (3V2a k v s - 2A 2 (v 2 + t>f )) + 2iV2a k \ 
l - ^2A 2 (v 2 + v 2 s ) - 3V2a k v s + yj (3y/2a k v s - 2A 2 (v 2 + v\ 



(43) 
(44) 



+ 2AV2a k \ 2 v 3 s (45) 



While in (36) we cancelled the singlet-doublet mixing in the CP-even sector by a particular 



choice of parameters, such a mixing still occurs in the CP-odd sector. 
Finally the charged Higgs mass matrix in the basis (h^, h^* = h^) is given by 



M 



/2 



i i 
i i i 



where = \g\v 2 . It contains one massless Goldstone mode, and one massive eigenstate 

M 2 H± = M 2 w+ 1 -X 2 (2v 2 s -v 2 ), 



1 

71 



IK + V) 



(46) 

(47) 
(48) 



The charged Higgs sector is identical to the MSSM one because S is electrically neutral. The 
matching can be easily done with the substitutions 



b -> 



1 

i „ 



(49) 



Finally, we must ensure that all physical square masses are positive, which is equivalent to 
checking that our solution is a minimum of the potential. Moreover the constraint on M^± 
implies that we are not breaking the U(l) em . Up to now we only prevented the origin to 
be a minimum solution. So such a requirement will impose further constraints on the free 
parameters, that can be summarized as follows: 

sign(a fc ) = -signOs), \a k \ < 2v / 2A 2 |v 5 |, (50) 

and one of the two following options 



a) v 2 s > \v\ 

b) \X 2 (v 2 - 2v 2 ) < M 2 



From now on we shall assume that the lightest CP even scalar is the one coupled to the W's, 
thus M 2 < M 2 . Moreover we want also M 2 > M§. This implies that the only available 
option is a). This fixes also the lightest CP-even Higgs couplings to gauge bosons to be exactly 
as in the SM. The lightest CP-even Higgs couples to the charged Higgs with 

2 A 2 v 2 

^hH+H 2cw ~YlM 2, 

z 

where in our notation cw = cos 9w, $w = sin 9w with 9w the Weinberg angle, and the coupling 
is normalized according to the conventions in [26]. The first term is the MSSM contribution 
for a = — 7r/4 = (5 — it/2, while the term in A is the NMSSM correction. 

We now consider the impact of relaxing the tan (3 = 1 condition in the present analysis. We 
assume as usual that the scalar singlet is very heavy, and kvg 3> Mz,cik, a \- Within such 
approximation it is easy to derive [28J that a ~ /3 — tt/2 still holds, and H is again essentially 
decoupled from WW or ZZ. On the other hand, one has 

Mi ~ Ml cos 2 (2/3) + -AVsin 2 (2/3) - ^-v 2 (A - k sin(2/3)) 2 . (52) 

2 fir 

Hence, the presence of the singlet can still give a negative contribution to the light Higgs boson 



mass. To prevent the latter negative contribution, we generalize eq. (31), and assume 



fc = A/sin(2/3) (53) 
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Figure 1: . Contour of = 125 GeV in the (tan/3, |A|) plane (solid line). The dashed line is 
the X-SUSY upper bound. 



so that 

Ml ~ Mf cos 2 (2/3) + -AV sin 2 (2/3) . (54) 

Here, we are interested in values of A and tan /3 that satisfy Mh ~ 125 GeV. In Fig. [l]we plot the 
corresponding two-dimentional curve of |A| versus tan/3. The black continuous line represents 
Mh = 125 GeV. In the A-SUSY theory [45J, A is increased so that the interaction becomes 
non-perturbative below the unification scale. However, A should not exceed ~ 2, otherwise 
non-perturbative physics would appear below 10 TeV, spoiling our understanding of precision 
electroweak data in the perturbative theory. The dashed line in Fig. [T] represents the A-SUSY 
upper bound. Then, only low values of tan /3 are allowed, in the range tan /3 < 8. In particular, 
tan/3 = 1 corresponds to the minimal A value (A ~ 0.72). 

Finally, the relevant quantities for the charged Higgs phenomenology in the general tan /3 case 
are 

\ 2 V 2 \ 2 V 2 

sin (2/3) 2 

A 2 v 2 

XhH+H- = cos(2/3) sin(a + /3) + 2c^ sin(/3 — a) — ?TT72 cos ( a + fl) sin(2/3) . (56) 



3.2 Neutralinos and charginos 

The soft SUSY breaking gaugino mass terms in the Lagrangian read 

£ = - 1 - (M1A1A1 + M 2 \\\\ + M3A3A3) . (57) 
In the basis ip° = (Ai, A|, h° d , h^, s), the resulting mass terms in the Lagrangian read 

£ = -^(/)W) + h.c (58) 
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where 



/ M 1 



V ... 



—MzSw cos P MzSw sin /3 
M 2 M z cw cos (3 —M z cwsm(3 


' 






-^§cos/3 



(59) 



/ 



Because of supersymmetry and electric charge conservation the chargino sector is the same as 
the MSSM up to the substitutions (49 ), in the gauge-eig enstate basis ^ = (W+, H+, W~, Hj) 
the chargino mass terms in the Lagrangian are 



-^chargino mass o ) ^-c^ + ^- c - 



(60) 



where, in 2 x 2 block form, 



with 



M, 



c 



X 





X 



_ M 2 
y/2M w cos /3 







y/2M w sin /3 



V2 



Xv c 



(61) 



(62) 



3.3 Radiative Higgs boson decays 

The model we have chosen to work with leaves the FP Higgs boson decays to WW* and ZZ* 
final states at tree level unaffected compared to the FP SM predictions. Although it is easy 
in this framework to decrease the coupling at tree level by choosing a different tan /3 and the 
Higgs mixing parameter a, in the FP Higgs scenario this will also suppress the induced 77 rate 
because the latter is dominated by the H^-boson loop. Because fermiophobia by itself is able 
to explain the observed deficit in WW* channel [T7] . our choice is motivated by a maximized 
77 rate. Therefore the deviations from the FP SM predictions may happen only due to extra 
particles in the loop. Because in the FP Higgs scenario the flavour physics constraints on 
charged Higgs masses are largely removed and chargino could be light, those particles can be 
as light as their present lower bounds from LEP II. 

The free parameters at the EW scale which are relevant for our analysis are the following: 
tan/3, the gaugino masses Mi and M 2 , the [i term given by [l = Xvg/V2, the sign(//Mx) and 
sign(/iM 2 ). In the present model, the mass of the charged Higgs is fixed once the value of \x 
is given, see Eq.fl47|). Moreover, we have chosen the convention of keeping M 2 positive, and 
allowing sign(/i) to vary. We have set M h = 125 GeV, M x = 100 GeV and y/s = 7 TeV. 

Then, tan/3, sign(/i) and M 2 are free parameters. We chose to re-express and M 2 as 
functions of two physical mass parameters: the charged Higgs mass (Mjj+) and the lightest 
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chargino mass (M +), as follows 



Mo 



M h + ( M h+ ~ M w) sin 2 (2/3) - M 2 cos 2 (2/3) 

71 



(63) 



M\ 



2 { M 2 + 



/J 2 ) + M^ ( 1 - cos 4/3) - 2/iM^ sin 2/3 



2(M 2 + - /x 2 ) 



(64) 



There are two different values of the gaugino mass M2, corresponding to the same lightest 
chargino mass. For convention M2 > and, for each sign(/x) and \n\ value, only one of the two 
solutions is allowed. Finally, we recall that the A parameter is determined by Mh and /3 as 



IAI 



M 2 h 



M 2 Z cos 2 (2/3) 



v sin(2/3) 



(65) 



In Table 1, we give some numerical values of the input parameters and the corresponding 
derived fundamental parameters, for M^ = 125 GeV and M H + = 400 GeV. 



Input parameters 


Derived parameters 


sign(//) 


tan (3 


M v+ (GeV) 


|A| 


H (GeV) 


M 2 (GeV) 


+ 


1 


200 


0.72 


290.8 


271.2 


+ 


5 


200 


1.40 


125. 


125.8 




1 


200 


0.72 


290.8 


186.8 




5 


200 


1.40 


125. 


151.3 


+ 


1 


400 


0.72 


290.8 


340.8 


+ 


5 


400 


1.40 


125. 


379.6 




1 


400 


0.72 


290.8 


390.6 




5 


400 


1.40 


125. 


383.9 



Table 1: Numerical values of the input parameters and the corresponding derived fundamental 
parameters for M h = 125 GeV and M H + = 400 GeV. 



The corresponding decay widths for the radiative induced decays h — >■ 77 and h — > Z*y of the 
lightest CP even Higgs boson h, in the framework of pure FP NMSSM model, are reported in 
Appendix. 

In order to avoid a large tree-level Higgs decay into an invisible sector [IB] , that would destroy 
the potential enhancement of the Higgs decay into 77 [18] . we will require that the mass of the 
lightest neutralino (M x o ), which is the lightest supersymmetric state in our scenario, is heavier 
than half of the Higgs mass. Then, due to i?-parity, all other Higgs decays into two generic 
neutralino states h —> XiXp including the case in which one or both are virtual states, will 
automatically vanish. In addition, we also require the lower bound on the chargino mass to be 
90 GeV. 
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Figure 2: Radiatively induced signal rates of 125 GeV FP NMSSM Higgs decays h — > 77 (left) 
and h — > Z7 (right) at the 7TeV LHC as functions of the lightest chargino xt mass (^ x + ) for 
several charged Higgs boson H + masses (M H + ) as indicated in figures and for tan/3 = 1. The 
SM central value prediction (dashed line) together with la uncertainty band A(SM) and the 
FP SM Higgs prediction are also shown. The lines above (below) the FP SM line correspond 
to sign(M 2 fi) > (sign(M 2 fi) < 0) for h — > 77 and to sign(M 2 fi) < (sign{M 2 jj) > 0) for 
h Z7. 



Taking into account the results obtained so far, we have computed the h — > 77 and h — > Z"f 
signal rates for the FP Higgs boson in the FP NMSSM. We present our results in Fig. [2] where 
we plot the signal rates of those processes as functions of the lightest chargino mass M + for 
different charged Higgs boson masses as indicated in the figure. The SM predictions together 
with their uncertainties and the FP SM predictions are also presented. The la (green) band 
corresponds to the theoretical uncertainty on the SM production cross section by gluon-gluon 
fusion. We have not included the uncertainty band on the SM FP line since the corresponding 
theoretical uncertainty due to the VBF cross section is quite small and can be neglected in this 
context. 

As in the MSSM case, the dominant contribution to SUSY contribution to h — > 77 and h — > Z7 
amplitudes comes from the charginos loop. From Figs. |2j [3j we can see that for fixed chargino 
mass there are always two solutions for the one-loop SUSY amplitudes corresponding to h — > 
77 and h — > Z7 decays. This can be understood as follows. For values of M v + below the 
interesection point with the FP(SM) line, the double solution is mainly due to the sign of \i 
that controls the relative sign of the SUSY amplitude with respect to the SM one. The lines 
above (below) the FP SM line correspond to sign{M 2 \i) > (sign(M 2 fi) < 0) for h — > 77 and 
to sign(M 2 fi) < (sign(M 2 fi) > 0) for h — > Zj. The dependence of the curves by sign(M 2 fi) 
can be understood from Eqs.( 66|67| ) in appendix. 
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Figure 3: The same as in Fig. but for a Mh = 125 GeV FP Higgs boson at the 8 TeV LHC. 



However, for values of M x + above the intersection point, the sign of /x is fixed and the double 
solution corresponds to the fact that in the loop run two non-degenerates values of heavy 
chargino states at fixed M +. These two values of heavy chargino mass, for a fixed light 



chargino mass, correspond to the two different solutions for the Mi parameter in eq. (64). 
The kink point corresponds to the case where the two solutions for the heavy chargino masses 
coincide. As we can see, there is a non-decoupling effect of the SUSY contribution to the loop 
77 and Zy decay amplitudes in correspondance of the maximum value for the lighest chargino 
mass. 



As seen in the Figs.|2j[3j the FP NMSSM signal rates can be both bigger or smaller than the FP 
SM predictions. For very light sparticles the total rate in 77 channel can even exceed the SM 
prediction. On the other hand, the present fits indicate that the LHC observes fewer 77 events 
than predicted by the pure FP SM [T5]. This result can be easily explained in the FP NMSSM 
since also rate reductions of as much as 50% are possible for the chosen parameters. The 
absence of points in the half-plane above (below) the FP(SM) line for the curve corresponding 
to M H + = 200 GeV in the case of h — >■ 77 (h — > Zy), is due to the lighest neutralino mass 
constraint M x o > M^/2 and depends on our choice for Mi = 100 GeV. 

At the 8 TeV LHC the predictions are qualitatively the same but numerically different. We 
present the rates for 8 TeV LHC in Fig. [3] for the same model parameters as in Fig. [2j 

If the top Yukawa coupling of the Higgs boson is not exactly vanishing, also gluon-gluon fusion 
production process will contribute to the Higgs production. In that case it is important to 
know our predictions for the FP Higgs branching fractions in our scenario. In Fig. [4] we plot 
the deviation of FP NMSSM Higgs boson branching fractions from the SM prediction for the 
previously specified parameters. The qualitative behaviour of branching fractions is the same 
as in previous figures, explaining our results. 
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Figure 4: Dependence of the FP NMSSM Higgs boson branching fraction over the SM value as a 
function of the lightest chargino mass M x + for Mh = 125 GeV, for different values of the charged 
Higgs boson mass M^j- and for tan (5 = 1. The dashed line corresponds to the FP SM scenario. 
The lines above (below) the FP SM line correspond to sign(M 2 fi) > (sign(M 2 fi) < 0) for 
h — > 77 and to sign{M 2 jj) < (sign{M 2 ji) > 0) for h — > Zj. 



Finally, in Figs. [5] and [6] we present the same plots as in Figs. [2] and |3J respectively, for 
tan/3 = 5. We set this value as an intermediate point between tan/3 = 1 and the A-SUSY 
upper bound tan/3 < 8 at M h = 125 GeV (see Fig. [I]). By increasing tan/3, the deviations from 
the SM FP Higgs predictions for production rates and BR's for h — y 77 and h — > Z7 are quite 
decreased. The largest effect is indeed achieved for tan/3 = 1. On the other hand, at tan/3 = 5, 
the largest SUSY contribution is obtained for a charged Higgs mass M H + ~ 400 GeV and a 
light chargino mass M x + < 150 GeV. The curves corresponding to M#+ = 200 are not present 
in Figs. [5] and [61 not being allowed for tan/3 = 5, because of the constraint M x o > M h /2. 



4 Discussion 



Apart from the Higgs boson phenomenology at the LHC discussed in the previous section, the 
FP NMSSM scenario has other important implications for SUSY phenomenology. As we have 
emphasized, the b — > S7 and B s — > fifi constraints on the charged Higgs mass are absent in this 
model. Thus the charged Higgs boson can be light and kinematically accessible at the LHC 
in the process pp — > H + H~ . The same may apply to other possible scalar and pseudo-scalar 
final states. While we have chosen to work with a particular model in which the Higgs boson 
is exactly SM-like, in general also other final states are possible, allowing to study the model 
parameters. However, because they are fermiophobic, their search strategy must be revised. 
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Figure 5: The same as in Fig. \^ but for tan /3 — 5. 



If neutralinos and charginos Xi are light? the dominant decay modes of all heavy Higgs bosons Hi 
could be into two Xi- m particular, the tree level decays of sr are induced by the A coupling. If 
the decay channels to sparticles are kinematically forbidden, the heaviest among A®, H, H ± 
will have tree level decays into the lighter ones and to (real or virtual) W's or Z's. Then the 
lightest of them, since it cannot decay into SM fermions because of fermiophobia, will have 
SUSY induced radiative decays giving in the final state the SM gauge bosons and fermions. 
We stress that, because of the decoupling induced by our values of a and (3, there are no 
trilinear vertices involving only one h and one of the scalars among A®, A\, H, H ± . Thus all the 
decays of heavier Higgs bosons are characterized either by large invisible branching fraction or 
multi-particle final states. Those decay signatures can easily be missed in present LHC searches 
explaining the absence of another Higgs-like resonance at higher masses. 

The second most relevant phenomenological implication of our framework concerns direct dark 
matter searches. In the MSSM the spin-independent dark matter scattering off nuclei is dom- 
inated by tree level Higgs boson exchange. In the FP Higgs case this process is suppressed. 
The dominant dark matter scattering process is through WW exchange at one loop level. This 
implies that the scattering cross section is suppressed by additional loop factor compared with 
the MSSM expectations. Scattering due to iy-loops is too weak [37J to be observed in the 
present stage of XENON100 



Arguably the biggest drawback of our scenario is the absence of explanation for the third gen- 
eration fermion masses. However, models of composite Higgs boson that explain naturalness 
of the electroweak symmetry breaking with new strong dynamics at 2-3 TeV scale do predict 
non-standard Higgs boson coupling to fermions [38J. Fermiophobia may be a feature of this 
framework. Supersymmetrizing the theory will stabilize the radiatively generated Yukawa cou- 
plings against new physics at high scales. While in non-SUSY case one expects large radiative 
corrections to Yukawa couplings proportional to log(M/ l /A), where A is the unknown scale 
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Figure 6: The same as in Fig. ^ but for tan /3 — 5. 



of new physics, in the SUSY version of FP Higgs those corrections will be at most of order 
log(M/i/MsusY), hence stabilizing the theory. 



5 Conclusions 

If there is a signal of a fermiophobic, or partially fermiophobic, Higgs boson with mass = 125 
GeV, the fundamental idea of supersymmetry, as it is implemented in the MSSM, is in trou- 
ble and must be revised. In particular, we have shown that the MSSM with vanishing or 
strongly suppressed Yukawa couplings is ruled out, independently of the particular supersym- 
metry breaking mechanism. Indeed, due to the absence of Yukawa couplings the usual (large) 
logarithmic corrections to the Higgs mass, induced by the scalar particles running in the loops, 
are absent and the upper bound on the Higgs mass is very close the Mz mass. 

In order to rescue supersymmetry, we show that a viable model beyond MSSM could be the 
NMSSM, where the absence of tree-level Yukawa couplings in the superpotential is guaranteed 
by the addition of a Z 3 discrete symmetry. The most relevant aspects of this scenario is that the 
SUSY naturalness criteria are automatically relaxed by a factor Ncyf/g 4 ~ 25, solving the little 
hierarchy problem and allowing sparticle masses to be naturally of order 2-3 TeV. Moreover, 
the usual flavor and CP problems are relaxed partly because of the absence of Yukawa couplings 
and partly for the possibility that the scalar partners are naturally heavy. 

In this framework, we consider the particular NMSSM case in which the mixing of the singlet 
with doublet Higgs fields is absent in the CP-even sector and at tree level the lightest Higgs 
boson is exactly SM-like. We analysed the predictions of this scenario for a Mh = 125 GeV Higgs 
at the LHC. We show that the predictions for the one-loop Higgs boson branching fractions 
and production rates in 77 and Z7 can be sizably modified with respect to the FP SM model, 
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allowing a better fit to present collider data. However, the tree-level Higgs decay channels into 
WW* and ZZ* remain unaffected if the mixing between the singlet and doublet Higgs fields is 
absent. Relaxing this last condition, and so adding a new free parameter, the Higgs coupling to 
weak gauge boson WW and ZZ can be modified, and a suppression of the rates for h — > WW* 
and h — > ZZ* with respect to the pure FP model expectations can be achieved. 

Finally, we would like to stress that the FP NMSSM offers a new arena for SUSY phenomenol- 
ogy at the LHC In particular, most of the previous analyses on SUSY particle searches should 
be revised in the light of the fact that the large top- Yukawa coupling is absent or strongly 
suppressed. In addition, the stringent constraints from Higgs mediated and other FCNC pro- 
cesses can be relaxed due to the absence or suppressed Yukawa couplings, allowing for a light 
charged Higgs boson phenomenology at the LHC. Moreover, the interplay between chiral and 
supersymmetry breaking suggests that if there is a new strong dynamics at the TeV scale, as for 
instance indicated by the large top-quark mass, this could also play a role in the supersymmetry 
breaking mechanism, opening the way to a new and exciting phenomenology at the LHC. 



Appendix 



Here we provide the analytical expressions for the one-loop radiative decays widths of h — >• 77 
and h — > Zj, where h is the lightest CP even Higgs boson, in the framework of pure FP NMSSM 
model. Following the results of Refs.[26], [H] we get 



r(h 77) 



+ 



a 2 G 



128\/27r 3 
2M W 



9hww^l{ T w) + 



E 



=1,2 



M, 



-9 h 



(66) 
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+ 
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64tt 4 



M^v H ±X hH+H - 
2c w M 2 
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9hwwAf(r w l , X w ) 
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2M W 

i=l,2\m=L,R ^xt 



E 



9hxtx-9™ xtx -A z 1/2 (r x l\ x +) 



(67) 



19 



where Gf is the Fermi constant, a the electromagnetic fine structure constant, and the nor- 
malized Higgs and Z couplings appearing above are given by [26] 



9h WW 



sin(/3 — a) 
1 



9 ktxj 



\[2s w 
1 



sin a VaU i2 + cos a V i2 Ui\ j . 
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cw V " 2 
1 / , I 



v i2 v il 



c w \ Sw 2 U " f/il) ' 



2r 2 
ZC W 



1 



H 



CW 



(6* 



where cw = cos 9w and sw = sin 6^, with 6 W the Weinberg angle, Xhh+h- is given in eq. (56), 
Uij and the matrix elements of the corresponding U, V matrices diagonalizing the chargino 



mass matrix X in Eq.(62) as WX.V 1 , and M Xi the corresponding eigenvalues. 



as n = Ml/ (AM 2 ), Xi = AM 2 /M 2 Z , with i = W,H + ,xt> while the functions A7 1/201) (a;) , and 



The other symbols appearing in the expressions of the widths in Eqs.(|66|) and (\67\j are defined 

as r, = M 2 /(AM 2 ), \ = 41 

A fi/2fi,i)( x >y) are g iven by 
• for h — > 77 

Aj /2 (x) = 2 [x + (x - 1) F(x) ] x~ 2 
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(69) 



for h — > 

A i/2( X >V) = h(x,y) - h{x,y) 



A$(x,y) 



h{x,y) 



c w U 3- 



' w 



-w 



h{x,y) + 



1 + 



'ir 



x l c 



where the functions Iip{x,y) are given by 
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(71) 
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with Fix) 



arcsm \/x 



for x < 1, and G(x) = y arcsin y^x for x < 1. The electromag- 
netic coupling constant a, appearing in eqs. (66) and (67), is evaluated at the scale q 2 
since the final state photons in the Higgs decays H — > 77 and H — > Z7 are on shell. 
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